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ABSTRACT 

The  Schwartz  Inequality  is  used  to  derive  the  Barankin  lower  bounds  on 
the  covariance  matrix  of  unbiased  estimates  of  a  vector  parameter.  The 
bound  is  applied  to  communications  and  radar  problems  in  which  the  unknown 
parameter  is  imbedded  in  a  signal  of  known  form  and  observed  in  the  pre¬ 
sence  of  additive  white  Gaussian  noise.  Within  this  context  it  is  shown  that 
the  Barankin  bound  reduces  to  the  Cramer-Rao  bound  when  the  signal  -to- 
noise  ratio  (SNR)  is  large.  However,  as  the  SNR  is  reduced  beyond  a  critical 
value  the  Barankin  bound  deviates  radically  from  the  Cramer-Rao  bound 
thereby  exhibiting  the  so-called  threshold  effect. 

A  particularly  interesting  signal,  which  has  been  widely  used  in  practice 
to  estimate  the  range  of  a  target,  is  the  linear  FM  waveform.  The  bounds 
were  applied  to  this  signal  and  within  the  resulting  class  of  bounds  it  was 
possible  to  select  one  which  led  to  a  closed  form  expression  for  the  lower 
bound  on  the  variance  of  the  range  estimate.  This  expression  clearly 
demonstrates  the  threshold  behaviour  one  must  expect  when  using  a  non¬ 
linear  modulation  system. 

Tighter  bounds  were  easily  obtained  but  these  had  to  be  evaluated  using 
numerical  techniques.  It  is  shown  that  the  side- lobe  structure  of  the  linear 
FM  compressed  pulse  leads  to  a  significant  increase  in  the  variance  of  the 
estimate.  For  a  practical  linear  FM  pulse  of  1  microsecond  duration  and 
40  megahertz  bandwidth  it  is  shown  that  the  radar  must  operate  at  an  SNR 
greater  than  10  dB  if  meaningful  range  estimates  are  to  be  obtained. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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I. 


INTRODUCTION 


An  important  problem  in  communications  and  radar  theory  is  the 
estimation  of  a  set  of  parameters  (0^,  0^,  ...»  ®m)  a  signal  which  has 
been  corrupted  by  additive  white  Gaussian  noise.  In  particular,  the 
received  waveform  is  assumed  to  be  of  the  form 

r(t)  =  s(t;9)  +  n(t)  |t|  =  T  (1-1) 

where  s(t;0J  is  a  known  function  of  t  for  each  value  of  the  m-row  vector 
0_.  In  the  radar  problem  0^  might  represent  an  unknown  time  delay 
(target  range),  an  unknown  doppler  shift  (target  velocity)  and  0^  an 
unknown  carrier  phase  angle.  The  noise  term  n(t)  represents  the 
Gaussian  white  noise  and  is  assumed  to  have  two-sided  spectral  density 
Nq  watts/Hz.  It  is  of  practical  and  theoretical  interest  to  determine  how 
well  a  given  estimation  scheme  can  perform  with  respect  to  estimating 
the  unknown  parameters.  In  this  respect,  Barankin  [l]  has  derived  a 
general  class  of  lower  bounds  on  the  moments  of  unbiased  estimators. 
Kieffer  [2]  has  used  the  Schwartz  Inequality  to  obtain  lower  bounds  on  the 
variance  of  an  unbiased  estimate  of  a  scalar- valued  parameter.  Applied 
to  the  pulse- position  modulation  communications  problem,  this  bound 
was  shown  to  yield  considerable  information  regarding  the  nonlinear 
modulation  threshold  effect  [3].  In  this  paper,  we  use  the  Schwartz 
Inequality  to  derive  lower  bounds  on  the  error  covariance  matrix  for 
unbiased  estimates  of  the  vector  parameter  0^  Then  we  specialize  these 
results  to  the  communications  and  radar  problem  as  formulated  in 
Equation  (1-1)  and  apply  the  bound  to  the  particular  problem  of  estimat¬ 
ing  the  range  of  a  target  using  a  linear  FM  waveform  over  an  incoherent 
channel.  It  is  shown  that  the  side- lobes  of  the  cor  responding  compressed 
pulse  significantly  affect  the  variance  of  the  estimate  of  the  time  delay. 
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H. 


THE  BARANKIN  BOUND 


Let  Q  be  a  sample  space  of  points  uo  and  let  P(u)/0)  be  a  family  of 
probability  measures  on  Q  indexed  by  the  parameter  9  taking  values  in  some 
index  set  it.  Assume  these  measures  have  a  density  function  with  respect 
to  some  measure  p,  i.  e.  ,  there  exists  a  function  p(aj  /0 )  such  that 

P(E/9)  =  ^p(w/0)dp(w) 


for  all  measurable  sets  E. 

Let  g(*  )  be  a  real  valued  function  defined  on  it  and  let  g(»  ) 
be  an  unbiased  estimator  of  g(9),  i.  e.  ,  §(•  )  is  a  real  valued,  measurable 
function  defined  on  Q  with  the  property  that 

Jg(uu)p(uu/0)dp(uo)  =  g(0)  (2-1) 

In  Appendix  A,  we  have  used  the  Schwartz  Inequality  to  show  that  the 

variance  of  the  estimator  g(.  )  when  0  is  the  true  value  of  the  unknown 

Z 

parameter,  denoted  o  (g),  is  bounded  below  according  to 

0 


{  2  a.tgO.)  -  g(0)]}2 


“  n 

r 

2  a  pfoi/gj) 

/ 

i=l  1 

/ 

p(d)  /0) 

p(uu  /0  )dp(uu ) 

which  is  valid  for  all  finite  families  (0.,a^).  It  follows  from  (2-2)  that, 
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(2-3) 
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a02(g) 


1.  u.  b. 

{  0  . ,  a .  } 

i  i  J 


{?  a.[g(e.)  -  g(0)]}  2 


n 

£  a  p(uu/e.) 

i=l  1 


p(ti)  /0) 


2 


p((JU  /0  )d(JL  (uu  ) 


where  the  1.  u.b.  is  to  be  taken  over  all  finite  families  of  9.e*n  and  real 

i 

a.. 

i 

Equation  (2-3)  is  the  Barankin  bound.  Barankin  has  shown  that 
this  is  the  best  possible  bound  in  the  sense  that,  for  each  0,  there  exists 
an  unbiased  estimator  that  achieves  it.  We  will  not  demonstrate  this  here. 

The  preceding  argument  has  assumed  no  structure  on  the 
parameter  space  tt.  We  now  assume  that  tt  is  the  m- dimensional  Euclidean 
space  E  of  vectors  9  and  we  will  use  (2-3)  to  derive  lower  bounds  for 

^  A  A  ^  A 

the  covariance  matrix,  E(9_  -  S_)  (9_-  9)  ,  of  unbiased  estimators  9_  of  0_. 
One  special  case  of  these  lower  bounds  will  be  shown  to  be  the  familiar 
Cramer- Rao  bound. 

Let  9  (u) )  be  an  unbiased  estimator  of  6_,  i.  e.  , 

E(9.)  =  J'Q.flu  )p(uu  /9)d|j(u))  =  £  , 

for  all  If  ^  is  an  arbitrary  m- vector,  then  g(uu  )  =  )  is  an 

unbiased  estimator  of  g(9^)  =_y/0^and  the  bound  given  by  Eq.  (2-2)  can  be 
applied  to  it.  This  yields  the  inequality, 


The  prime  denotes  matrix  transpose 
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2.„, 

CTe 


=  }r'E(0_  -  _0  )  (©  -_0)'  y_  r 


n  , 

{y '  s  ai(ii  -  !), 
i=l 


p(w  /0  )d(j(uo ) 


(2-4) 

and  this  is  valid  for  all  finite  sets  of  vectors  9 .  in  E  and  all  real  a..  For 

— i  m  i 

a  given  set  of  0  .  and  any  set  of  a.'s  we  choose,  (2-4)  yields  a  lower  bound 
ona02(g). 

It  is  therefore  reasonable  to  seek  that  set  of  a.  which  leads  to  the 

i 

least  upper  bound  of  (2-4).  The  maximization  is  performed  in  detail  in 
Appendix  B  where  it  is  shown  that  the  optimum  set  of  a^  for  a  fixed  set  of 
JK  yields  the  bound 

E(9  -0)(0-0)/  =  A"1  +  (§-  A_1A)  A"V  -  A_1A)7  (2-5) 


where 


A  =  B  -  A'A_1A 


(2-6) 


A  •  P(m/e)d|i(u))  (2-7) 

J  .  i  3 


A  = 


-I' 


i,  j  =  1,  2,  .  .  .  m 
?  to  pfa  /e  j - .  p(ai/e)dM(uj )  (2-8) 


p(uo/e_) 


i  -  1,  2,  iii , 


k  =  1,  2,  .  .  .  ,  n 


A  matrix  A  is  said  to  be  less  than  or  equal  to  a  matrix  B,  A  <B,  if 
B  -  A  is  non- negative  definite. 
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B  = 


p(uj  /9)dp(uu ) 


(2-9) 


/ 


P(»J  /  Q_i ) 

P(u>  W) 


p(m  /9k) 

P(ju  /9 ) 

i,  k  =  1 ,  2,  .  .  ,  n 


§  =  [9_1,  £2,  .  .  .  ,  e_n] 


(2-10) 


Notice  that  0  i  =  1,  2,  .  .  .  ,  m  refers  to  the  i  component  of  the  true 

parameter  value  0_,  while  0^,  i  =  1,2,  n  refers  to  an  arbitrary  value 

of  the  parameter  0_  other  than  its  true  value.  These  vectors  can  be  chosen 

at  will  and  the  number  of  them,  n  is  arbitrary.  For  reasons  which  will 

n 

become  clearer  in  the  sequel,  we  refer  to  as  the  set  of  test  points. 

If  we  pick  n  =  0,  the  second  term  on  the  right-hand  side  of 
Equation  (2-5)  vanishes  and  we  obtain  the  Cramer- Rao  bound. 


E(0_-  0_)(0_  -_0)'a  A-1 


For  n  >  0,  we  obtain  other  bounds  that  are,  in  general,  an  improvement 
on  the  Cramer- Rao  bound.  This  follows  from  the  fact  that  A  is  positive 
definite  (Appendix  B)  which,  in  turn,  implies  that  the  matrix  ($  -  A  A) 
A~*(§  -  A”^A)/  is  at  least  positive  semi- definite. 
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in. 


APPLICATION  TO  COMMUNICATIONS  AND  RADAR 


It  is  of  interest  to  specialize  the  results  just  obtained  to  the 
problem  of  estimating  the  parameters  of  a  signal  corrupted  by  additive, 
white  Gaussian  noi  se.  More  precisely,  we  are  given  a  stochastic 
process  r(t)  of  the  form, 

r(t)  =  s(t,9j  +  n(t),  \t\  =  T  (3-1) 

where  s(t,  0_)  is  a  known  function  of  t  for  each  6_  and  the  n(t)  is  zero  mean 
white  Gaussian  noise  process  with  the  covariance  function  E[n(t)n(s/)]  = 
NQ6(t-s/).  In  addition  we  are  given  an  unbiased  estimate  of  9_,  i.  e.  ,  a 
measurable  vector  valued  function,  0_(«  ),  defined  on  the  space  of  all 
sample  functions  r(»  )  and  satisfying 

E(|)  =  0_ 

We  want  to  obtain  a  lower  bound  for  E(9^  -  0_)(9_  -  G)7.  In  this  case  Q  is  the 
space  of  all  functions  r(-  )  defined  on  -  T  £  t  £  T  and  p(uj  /Q)  is  the  probability 
density  of  r(*  )  relative  to  the  measure  defined  by  the  white  noise  n(t)  alone, 
i.  e.  , 


p[r  (•  )/0j  =  exp  i —  J  r(t)s(t,0_)  dt  -  —  J  sZ(t,  0_)dt^ 

O  O  rp 


In  order  to  calculate  the  A,  B,  and  A  matrices,  it  is  sufficient  to 
evaluate  the  function 


GQ'.e^O)  =  f  P^/aXp/fi.*)  p(u,/0)d|a(u,) 

u  p  («i/e_) 
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For  the  problem  at  hand,  this  function  is  given  by  the  expression 


1  T  1  T 
exp  — J  r(t)s(t;0_' )dt  exp  — J  r(t)s(t;0_'/ )  dt 

g(9  ',0";9)  =  E  '  °  _T  °  "T 

exp  [  r  (t)s (t;9_)dt 

o  -T 

T  T 

exp  TJ—  J  s2  (tjG^7)  dt  exp  275-  J  s2(t;0^)dt 


o  -T 


o  -T 


,  i 

exp  J  _  ss  (t;£)  dt 


=  E 


o  -T 


T  T 

exp  N~  J*  n(t)  s (t;9_)Jdt  exp  J  n(t)js (t;£") 

-s(t;9)]dtj 


T  T  *i  / 

exp J  js(t;9_)-s(t;0  ')J2  dt  exp  J  [s(t;0_)- s(t;0_'  )J2  dt 


T 
o  -T 


Now  it  is  easy  to  verify  that 


where  x  is  a  zero  mean  Gaussian  random  variable  with  variance  a 
our  case 


2 

x 


In 
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X 


1 


T  .  T 

f  n(t)[s(t,£')  -  s(t,0_)]  dt  +  -rj-  J  n(t)[s(t,_e" )  -  s (t,  0_)J  dt 
o  -T  o  -T 


2 


a 

X 


T 

J  [s(t,2/)  +  s(t,  0_")  -  2s(t, 


dt 


Theref  ore, 


T 

G(Q_' ,  Q_"  ;QJ  =  exp  J  [s(t,e_')  -  s  (t,  9_)][s  (t,_©  "  ) 

o  -  T 


s(t,_0  )]  dt  J- 


(3-2) 


Using  Eq.  (3-2),  it  is  possible  to  calculate  explicit  expressi  ons  for  A,  B, 
and  A  as  follows, 


J _  p*  d  s(t,9.)  5  s(t,&)  ,, 

No  -  t  59.  ae. 

i  k 

i,  k  =  1,  .  .  .  ,  m 


(3-3) 


A  = 


J  ds(t,  6_) 

J  59. 

-T  i 


[stt.e^)  -  s(t,e_)]  dt 


k  =  1,  2,  ....  n 
i  =  1,  2,  .  .  .  ,  m 


(3-4) 
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B  =  exp 


(3-5) 


1  ^ 

-j-p  J[s(t,0_.)  -  s(t,  0_)][s(t,  e_k)  -  s(t,9_)]  dt 
o  -  T 

i,  k  =  1 ,  2,  .  .  ,  n 

Using  these  relations  the  Barankin  Bounds  applied  to  the  communications 
and  radar  problem  becomes 


E(0-9  )(9-9  ) '  ^  A'1  +  ($-A'1A)A"1($-A'1A)/ 


(3-6) 


where 


A  =  B  -  A  'A_1A 


$  ~  [ii*  '  •  *  ’in] 


(3-7) 

(3-8) 
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IV. 


LINEAR  FM  RANGE  ESTIMATION  ACCURACY 


In  the  radar  ranging  problem  a  known  signal 


a(t)cos(uu  t  +  cp (t ) )  |t|  =  T/2  (4-1) 


is  reflected  from  a  stationary  point  target,  distorted  in  phase  and 
observed  in  the  presence  of  noise  at  the  receiver.  The  received  wave¬ 
form  can  be  written  as 


a(t-  t)cos[o)  t  +  cd ( t  —  t)  +  y]  +n(t)  \  t|  =  T  (4-2) 


where  the  parameter  t  represents  the  time  delay  or  radar  range  while 
Y  represents  a  random  phase  shift  (uniform  in  [0,  2tt]).  A  typical  radar 
receiver  attempts  to  estimate  the  unknown  time  delay  by  using  a  likelihood 
signal  processor.  Since  the  radar  must  operate  in  a  variety  of  noisy 
environments,  (Nq  large  and  small),  it  is  important  to  know  how  well 
a  given  signal  is  likely  to  perform  for  a  range  of  signal-to- noise  ratios 
(SNR).  By  applying  the  Barankin  bounds  to  the  above  signaling  problem, 
we  can  determine  the  lowest  values  of  SNR  above  which  the  radar  must 
operate  to  give  acceptable  performances.  These  remarks  will  become 
clearer  as  the  analysis  proceeds.  From  Eq.  (4-2)  we  see  that  for  the 
radar  ranging  problem  the  parameter  vector  is 


(4-3) 


The  linear  FM  waveform  has  been  the  subject  of  considerable 


interest  in  both  practice  and  theory.  In  the  remainder  of  the  paper  we  shall 
assume  that  the  signal  has  the  linear  FM  phase  function  so  that 
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(4-4) 


cp(t)  =  j- 


t 


2 


where  U  =  2ttWT,  W  is  the  signal  bandwidth  and  T  the  time  duration.  In 
addition  we  shall  assume  that  the  amplitude  modulation  a(t)  rises  very 
rapidly,  but  smoothly,  to  the  constant  value  ^2E/T  when  |t|  =  T/2 
and  that  it  is  zero  when  |t|  >  T/2.  Figure  1  illustrates  a  typical  a(t). 

Then  the  signal  we  are  dealing  with  has  the  functional  form 

s(t;9 )  =  a(t-T)cos[u)  t+y(t-T)  +y]  (4-5) 

where  u)  represents  the  carrier  frequency  in  radians/sec. 

In  Appendix  C  we  have  carried  out  the  detailed  manipulations  which 
lead  to  the  Barankin  bounds  for  the  linear  FM  signal.  It  is  shown  that  the 
variance  of  the  normalized  delay  error  is  lower  bounded  by  the  following 
expression 


E(t  -t 

o 


E/N 


12  I  /A-l 

— 2~  +  cp  A 


(4-6) 


where 


(4- 7a) 


(4- 7b) 
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Fig.  1.  A  typical  "on-off"  amplitude  modulation. 
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1 


cos  (AYi)sin  [^AT.fl  -1At.|)] 


b..  =  1  - 


U 

2 


At. 

1 


- - cos(AY.)sin[^  At  (1  -  |  At  |  )]  (4-  7c) 

Hat.  J  *  J  J 

2  J 


+  — -  co  s  (Ay-- Ay. )  sin  [H  (AT  At.  )(1  -  |  At  .- At.  1  )] 

H  (At. -At.)  1  •*  Z  1  J  1  J1 

2  1  J 


Hi 


and 


are  vectors  whose 


.th 

j  row  is 


aij  =  *±2  ‘  a?- *cos  (AYj)cos  (§■  ATj{1'  lATjl^ 

(4-  7d) 

+  — ;  -1-? - cos(AY-)  sin[H  At  (1  -  |at  | )] 

H  At.  J  *  J  J 

2  J 


where  the  ±  signs  correspond  to  At  >  0 


a 


1 


(±2  - 


1 


At. 

J 


)cos[H  AT(1  -  |at.|)] 
c  J  J 


4-7p.) 


AT.  =  (1.  -  To)/T 


(4- 7  f ) 
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Notice  that  the  bounds  depend  only  on  the  delay  and  phase  differences 

T.  -  T  ,  Y  .  -  Y  and  not  on  the  true  values  of  delay  and  phase  T  ,  Y  . 
JOJO  <  <  oo 

Since  we  know  that-T/2  =  t  =  T/2,  then  we  need  study  only  the  range 
-  T  =  T-tq  =  T,  or  equivalently,  -  1  =  At^  =  1.  It  is  reasonable  to  assume 
that  the  phase  is  uniformly  distributed  in  [0,  2ir]  so  that  the  test  points, 
AYj,  can  be  chosen  equally  spaced  in  the  [0,  2tt]  interval.  We  are  free  to 
put  the  test  points  anywhere  we  like  within  the  above  intervals,  and  for 
each  set  of  these  points  we  get  a  lower  bound  on  the  variance  of  the 
normalized  delay  error.  To  gain  some  insight  into  how  we  might  locate 
these  points,  let  us  first  examine  the  case  in  which  there  are  no  phase 
test  points.  As  a  result  the  following  simplified  equations  are  obtained 


a 


[^At  (1  -IAt  1)] 

J  J 


(4- 8a) 


^  j 


(4- 8b) 


where  the  ±  sign  corresponds  to  At.  <  0. 


(4-9) 


+ 


— - - -  sin[£  (at.  -  AT.)(1  -  j  At±  -  At.  \ )] 

(At.  -  At.)  j  j 

c  1  J 
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It  appears  that  the  significant  terms  depend  on  the  quantity  sin[~  x(l 

2“X  which  is  well-known  to  be  the  response  of  the  filter  which  is  matched  to 


the  linear  FM  waveform.  This  filter  response  is  used  in  the  implementa¬ 
tion  of  the  maximum  likelihood  receiver.  It  is  interesting  that  the  Barankin 
bounds  should  lead  to  terms  which  depend  on  this  receiver  structure.  We 
know  from  physical  arguments,  that  the  likelihood  processor  deviates  from 
the  Cramer-Rao  performance  when  the  noise  becomes  strong  enough  to  cause 
the  peak  detector  to  lock  on  one  of  the  sidelobe  peaks  of  the  matched  filter 
output.  Since  the  Barankin  bounds  appear  to  be  taking  this  sidelobe  behavior 
into  account,  we  shall  use  this  intuitive  knowledge  to  justify  locating  the 
delay  test-points  at  the  sidelobes  of  the  sin^-x(l  -|x|)/^-x  function. 


As  an  example  let  us  choose  two  test  points 


ATj  =  -  At2  =  6 


where  6  is  the  location  of  the  first  positive  sidelobe.  For  large  time- 
bandwidth  products  we  can  approximate 


and  therefore  we  choose  ^  6  =  5*11/2  which  leads  to  6  =  2.  5/WT.  Therefore 
the  approximation  is  fairly  good  for  WT  ^  50.  Then 


cos  — -  6(1  -  6)  «  cos  ^6=0 


and  Eqs.  (4- 8a),  (4-9)  reduce  to 
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an  =  +  WT/6.  25ir 

al2  =  '  WT/6-  25ir 
a2l  =0 
a22  =  0 


bil  =b22  =2 


1  - 


— —  sin  -^6(1  -  6  ), 

U  c  2 

2  6 


=  1. 746 


U  , 
2  6 


5tr 

2 


bl2  b2 1 


1  -  — sin  £  6(1  -  6) - - -  sin  £  (-6)(l-6) 


*<-*> 


+  -jjg-  sin  u6(l  -26) 


=  0.746 


- 


5ir 

2 


Substituting  these  values  in  Eqs.  (4-7a)  and  (4- 7b)  we  find  that 


12  WT 

ep,  =  -cp?  =  6  -  - - T  =  2.485/WT 


(2uWT) 


'll  22 


A77  =  exp  1.746 


12  /  WT 

!  1 


o  1N  o  (2uWT)2  y6*25Tr 


=  exp  1.  746  -  0.  00078 

o  o 


exp  1.  746 

o 
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A12  =  A21  =exP  °'746  W  +  #" 

O  Q 


12  _  /  WT 

2  1  6. 25u 


(2irWT) 


=  exp  0.  746  -£-  +  0.  00078  §- 
N  N 

o  o 


;  exp  0.  746 


N 


In  this  case  the  A- matrix  is  easily  inverted  to  give 


-1 


"Z  ~T 


A11  ‘  A12 


A11  "  A 12  \'A12  A11 


and  the  quadratic  form  appearing  in  Eq.  (4-6)  reduces  to 

T  - 1  2 

cp  A  —  2  cp  ^  /  (A  ^  j  ~  A  ^  ^ ) 


=  2 


2.485 

WT 


exp  1.  746 -  exp  0.  746 

o  o 


12.38  exp  (-1.746^=-  ) 

o 

1  -  exp  (-  — ) 


17 


251 W 


Substituting  this  result  into  Eq.  (4-6)  we  obtain  a  lower  bound  on  the  variance 
of  the  normalized  delay  error,  namely 


1  it/"  ^  > 

IT  E(t'to)  = 


(WT)2 


N 


N 


where 


40.  6  „ E  exp  (-1.746  E 


N 


N 


N 


=  1  + 


1  "  exP 


When  ^  is  large,  T  (  ^  )  =»  1  and 
o  o 


E  (t-to) 


2  > 


4 —  (WT)2 


1 

E 


N 


(4- 10a) 


(4. 10b) 


which  is  the  Cramer-Rao  lower  bound  for  the  linear  FM  waveform.  It  has 

been  shown  [4]  that  the  performance  of  the  matched  filter  receiver  achieves 

this  bound  when  -rj -  is  large,  but  it  is  not  known  when  the  performance 

begins  to  deviate  from  it.  Letting^ -  =*  0  in  Equation  (4- 10b)  we  see  that 

E  ^o 

T  (~ — )  =*»  and  the  so-called  threshold  effect  has  taken  place.  For  this  reason 
tN  ©  e  E 

we  refer  to  T  (rrz - )  as  the  threshold  function.  The  value  of  — - at  which 

N0  N0 

the  tighter  bound  of  Equation  (4- 10a)  deviates  from  the  Cramer-Rao  bound 
is  referred  to  as  the  threshold  operating  point.  Therefore,  Eq.  (4- 10a)  gives 
the  first  analytical  expression  which  determines  the  best  performance  the 
matched  filter  can  ever  achieve.  Expressed  in  db  the  lower  bound  on  the 
normalized  delay  variance  is 
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101°S10  {^2 •  Eft-r/j  = 


10  log  10 


(2ttWT)2 

12 


_E 

N 


o 


+  10  log 


(4-11) 


In  Figure  2  we  have  plotted  the  bound  for  the  case  WT  =  100  and  note  that 
changing  the  time-bandwidth  product  merely  shifts  the  curve  by  a  constant. 

The  figure  clearly  shows  the  threshold  effect  which  we  referred  to. 

One  must  be  careful  to  give  the  correct  interpretation  to  the 
bound  described  by  Eqs.  (4- 10b)  and  (4-11).  It  states  only  that  the  perform¬ 
ance  of  any  radar  which  uses  a  linear  FM  waveform  to  estimate  range  can 
be  no  better  than  that  specified  by  Eq.  (4-11).  For  example  if  WT  =  100 
and  the  input  SNR  is  2.  5  dB,  then  the  normalized  range  variance  must  be 
something  greater  than  -41.  0  dB  which  is  7  dB  larger  than  that  predicted 
by  the  Cramer-Rao  bound. 

The  bound  can  also  be  used  to  give  a  lower  bound  on  the  input 
SNR  which  must  be  used  to  give  acceptable  performance.  The  Barankin 
bound  deviates  from  the  Cramer-Rao  bound  when  the  noise  at  the  sidelobe 
peaks  becomes  significant.  This  means  that  radar  will  be  making  large 
errors  in  range  and  the  performance  becomes  unacceptable.  If  we  locate 
the  threshold  operating  point  at  the  value  of  the  input  signal- to- noise  ratio 
at  which  the  Barankin  bound  deviates  from  the  Cramer-Rao  bound  by  1/2  dB, 
then  if  only  the  two  largest  peaks  are  considered,  Fig.  2  shows  that  the 
radar  must  be  operated  at  an  input  SNR  greater  than  6  dB.  This  is  not  a 
useful  result  from  a  practical  point  of  view  since  the  radar  usually  operates 
above  13  dB,  the  signal- to-noise  ratio  at  which  target  detection  occurs  with 
a  suitably  small  false  alarm  probability.  Since  the  curve  gives  only  an 
upper  bound  on  performance  we  cannot  say  conclusively  that  the  practical  radar 
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Fig.  2.  Lower  bound  on  range  variance  using  two 
side -lobe  test  points. 
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operating  at  13  dB  input  SNR  will  always  avoid  the  threshold  effect.  However, 
we  can  obtain  a  tighter  set  of  lower  bounds  by  increasing  the  nurrber  of  delay 
and  phase  test  points.  In  these  cases,  the  A  matrix  cannot  easily  be  inverted 
analytically  and  we  must  resort  to  numerical  techniques. 

In  Figure  3  we  show  the  effect  of  taking  more  of  the  sidelobe  peaks  into 

account  for  a  waveform  having  a  time -bandwidth  product  WT  =  40.  We  have 
plotted  the  lower  bound  for  the  cases  where  the  test  points  are  located  at  the 
first  two  negative  peaks,  then  consider  the  additional  next  two  positive  side- 
lobes,  etc.  ,  until  all  of  the  peaks  are  accounted  for.  When  all  the  sidelobes, 
are  considered  a  considerably  tighter  bound  is  achieved. 


In  the  next  experiment  we  located  the  delay  test  points  at  the  maximum 


number  of  sidelobe  peaks  and  selected  the  phase  test  points  equally  spaced  in 


the  [0,  2tt)  interval.  That  is,  if  M  is  the  number  of  phase  test  points,  we  set 


2t t 


i  =  1,  2,  .  .  ,  (M+l).  Figure  4  shows  a  sequence  of  curves 


AYi  =  (i-D 


(M+l) 


which  demonstrates  the  effect  of  increasing  the  number  of  phase  test  points. 
Numerical  results  were  obtained  for  several  time -bandwidth  products  up  to 
WT  =  40.  It  was  found  in  all  cases  that  no  significant  change  could  be  ob¬ 
served  when  more  than  four  phase  test  points  were  used,  and  that  the  lack 
of  phase  knowledge  results  in  a  loss  in  potential  performance  of  about  5  dB 
in  all  cases. 

In  a  practical  radar  system  evaluation,  one  would  like  to  relate  these 
bounds  to  potential  range  accuracy.  The  range  of  a  target,  R,  is  related  to 
the  round  trip  time  delay  as  r  -  — c  T— 

where  c  is  the  velocity  of  light  (3  x  108  meters/sec)  and  therefore,  the  stand¬ 
ard  deviation  of  the  range  estimate,  <JR,  is  given  by 


(4-12) 
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NORMALIZED  DELAY  ERROR  (dB) 


Fig.  3.  Lower  bound  on  range  variance  using  additional 
side-lobe  test  points  and  no  phase  test  points. 
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NORMALIZED  DELAY  ERROR  (dB) 


13-42-11943(1)1 


-40 


-30 


-20 


-10  — 


0  — 
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CRAMER -RAO 
LOWER  BOUND 


M  =  0 


M  =  NUMBER  OF  PHASE 
TEST  POINTS 
N  =  NUMBER  OF  DELAY 
TEST  POINTS 


40,  N  =  40 


L 


4  6  8 

INPUT  SNR  E/N„  (dB) 


10 


12 


14 


Fig.  4.  Lower  bound  on  range  variance  using  all  the 
side -lobe  test  points  and  additional  phase  test  points. 
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Experiments  have  recently  been  conducted  which  use  a  linear  FM  waveform 
of  40  microseconds  duration  and  1  megahertz  bandwidth.  Therefore  the 
Op  corresponding  to  this  pulse  can  be  computed  using  Equation  (4-12) 
and  the  data  used  in  plotting  Figure  4.  The  results  are  shown  in  Figure  5. 
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Fig.  5.  Lower  bound  on  range  error  for  a  1  megahertz 
40  microsecond  linear  FM  pulse. 
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V. 


CONCLUSIONS 


A  simpler  derivation  of  the  Barankin  bound  for  the  variance  of 
unbiased  estimation  of  vector  parameters  has  been  derived.  Under  the 
assumption  that  the  parameters  are  estimated  on  the  basis  of  a  signal 
imbedded  in  additive  white  Gaussian  noise,  we  have  reduced  the  bound 
to  a  form  which  is  suitable  for  numerical  evaluation.  It  is  shown  that  the 
variance  depends  on  two  terms;  the  first,  which  is  the  Cramer-Rao  bound, 
involves  the  structure  of  the  signal  in  the  vicinity  of  the  true  parameter 
value,  (fine  structure),  while  the  second  term  takes  into  account  the  system 
performance  at  parameter  values  far  removed  from  the  true  values  (gross 
structure).  It  is  this  remarkable  property  of  the  Barankin  bounds  which  is 
so  useful  for  evaluating  system  performance.  Applied  to  the  problem  of 
estimating  target  range  using  a  linear  FM  waveform,  we  see  that  the  fine 
structure  determines  the  performance  in  the  large  SNR  region  while  the 
gross  signal  structure  dominates  when  the  SNR  is  small.  The  bounds  are 
extremely  useful  from  a  systems  design  point  of  view  since  they  yield  the 
SNR  at  which  the  transition  in  performance  takes  place.  Usually  the 
performance  is  acceptable  when  the  fine  structure  is  predominant.  For  the 
waveform  under  study,  it  was  shown  that  the  radar  should  operate  above 
10  dB.  Below  this  value,  degradation  in  performance  was  rapid.  Although 
adequate  detection  capability  requires  an  input  SNR  above  13  dB.  ,  there  is 
no  reason  why  the  system  could  not  shift  to  lower  SNR's  when  the  target  is 
being  tracked.  The  bound  could  then  be  used  to  give  lower  bounds  on  this 
tracking  SNR. 

It  is  interesting  that  the  bounds  derived  here  can  be  related  to  the 
side- lobes  of  the  linear  FM  matched  filter  output.  It  has  long  been  a  design 
criterion  to  use  waveforms  which  have  a  low  side- lobe  structure.  The 
Barankin  bounds  lead  to  the  first  analytical  evidence  that  this  is  indeed  a 
good  design  criterion. 
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Conceptually  the  effect  of  an  unknown  doppler  shift  can  easily 
be  incorporated  into  the  analysis  although  the  mathematical  manipulations 
become  somewhat  more  complicated.  However,  it  would  be  of  interest  to 
see  how  much  the  lack  of  knowledge  of  the  target's  velocity  costs  in 
relation  to  the  accuracy  of  the  range  estimate. 

Currently  under  investigation  is  the  derivation  of  an  upper  bound 
for  the  variance  of  the  delay  error  when  the  range  is  esimated  using  a 
matched  filter  receiver.  In  conjunction  with  the  lower  bounds  presented  here 
the  tools  for  the  thorough  analysis  of  the  ranging  performance  of  a  radar 
will  then  be  availab  le. 
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APPENDIX  A 


The  following  is  an  informal  derivation  of  the  Bararikin  bound. 

Let  0  be  a  sample  space  of  points  id  and  let  P(u) /0 )  be  a  family  of  probability 
measures  on  Q  indexed  on  the  parameter  0  taking  values  in  some  index  set 
ir.  Assume  these  measures  have  a  density  function  with  respect  to  some 
measure  |a,  i.  e.  ,  there  exists  a  function  p(uu  /0  )  such  that 

P(E/0)  =  ^p  (id  /0  )d|_i  (uj  ) 


for  all  measurable  sets  E. 

Let  g(*  )  be  a  real  valued  function  defined  on  ir  and  let  g(«  )  be  an 
unbiased  estimator  of  g(0),  i.  e.  ,  g(*  )  is  a  real  valued,  measurable  function 
defined  on  fi  with  the  property  that 


t^g(w  )p(u> /0)dn(u) )  =  g(0)  (A-  1 ) 


We  now  obtain  a  lower  bound  on  the  variance  of  any  such  g(.  ). 

For  any  finite  set  of  points  0  .eir,  i  =  1,  .  .  .  ,  n  and  any  set  of  real 
numbers  a^,  i  =  1,  .  .  .  ,  n,  it  follows  from  (A-l)  that 


n  N 

t"g(uu)E  a.p(u)  /0i)dta((D)  =  £  aig(ei)  (A-2) 
i  =  1  i  =  1 

N 

We  subtract  a.g(e)  from  both  sides  of  (A-2)  to  obtain 


n 

J[g(u)  )  -  g(G  )]  £  a^p((D  /0i)d|i  (u) ) 

i  =  1 


[g  (e±)  -  g(e )] 


(A- 3) 
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n  and  all 


an  equation  valid  for  fKeir,  i  =  1,  .  .  .  ,  n,  all  real  a^,  i  =  1,  .  .  .  , 

0  e  tt. 

Equation  (A- 3)  can  be  rewritten  in  the  form. 


n  . 

E  a^p(uu/0.) 

J[g(w)  -  g(B  )]  - yQ) — p(w/0)d4(u)) 

2 

n 

s  a  -  [g  (0  - )  -  g(e )] 

i  =  1 

L  i 

(A- 4 

Applying  the  Schwartz  inequality  to  the  left-hand  side  of  Eq.  (A-4),  we 
obtain 


y  fg(0i)-g(9  )1(^  ]  ^[g(w  )  -g(e)]2  p(w  /e  )du(«j) 

i  =  1  L  I  | 


n 

2 

[ 

E  a  p(u)/0.) 

•  _  I  1  1 

J 

1  -  I 

p(io  /0  )d|a(ti) ) 

J 

p((l)  /0) 

i 

The  first  term  on  the  right-hand  side  of  Eq.  (A- 5)  is  the  variance  of  the 

2 

estimator  §(•  )  when  9  is  the  true  value  of  the  unknown  parameter,  aQ  (g). 

0 

We  now  have  the  inequality, 


|  £  ^a.[g(0.)  -  g (9 )]  | 


n 

E  a.pfuu/©.) 
i  =  1 


p(ti)  / 9 ) 


2 


p(uu  /0)d|j(u) ) 


(A- 6) 


valid  for  all  finite  families  {0^,  a^}. 
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APPENDIX  B 


In  Section  II,  we  showed  that 


n 


{y'  £  a.(0.  -  9)} 

y{E(@  -  e)(0-  e)'}*  =  r-n  1  =  1 - 


/ 


S  p(to  /e^ 

1-1  p(uu  /0) 


p(u)  /0)d|a((D) 


This  lower  bound  is  valid  for  every  set  of  0.eE  and  all  real  a..  We  now 

1  —l  m  i 

assume  a  fixed  set  of  0^  has  been  specified  and  then  seek  the  set  of  a^  which 
leads  to  the  least  upper  bound  of  (B-l).  To  do  this  we  now  specialize  some 
of  the  0,'s  anda.'s,  As sume  n  >  m  and  define 

—l  l 


0 .  =0  +  e.e. 
— i  —  l—i 


— m+1 


=  0 


i  =  1,  .  .  .  ,  m 


where  e.  /  0  is  a  real  number  and  e.  is  the  1 
1  —l 


.th 


unit  vector  in  E 


m 


i.  e. 


♦ 


.th 


place 


31 


Further  define 


X. 


m  X . 

m+1  i=l  e. 

1 


for  arbitrary  real  \^'s. 

The  remaining  6_^'s  and  a^'s,  i=m+2,  .  .  .  ,  n  are  left  arbitrary. 
With  these  definitions,  we  can  write 


!n  \  2  n 

T,  a  pfuu/e^l  =  E  aiakp((ju/0_i)p(au/e_k) 
i=l  )  i,  k=l 


m+1 

=  £  aiakp(au  /oi)p(«j/e4c) 

i,  k=l 


n 

=  E  a  a  p(cu /e  )p(uo /0,  )  + 

i,  k=m+2  1  K  _1  -K 


m+1  n 

+  2  E  E  a  a  p(uu /0  )p(w/9,  ) 

i=l  k=m+2  1  *  -1 
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m 


=  E 
i,  k=l 


p(a)  /9_  +  e^)  -  p(o)  /§_) 


0 . 

i 


p(uo  /9  +  Gj^)  -  p(u)  /e) 


m 
+  2  E 
i=l 


n 

E 

k=m+2 


Vk 


p(co/0_)  +  e^.)  -  p(uu/0_) 


e . 

1 


p(o)  /e_k) 


n 

+  s  a  a,  p(u)/0.)p(uo/e,  ) 

i,  k=m+2  1  k  -1 


It  now  follows  that  the  denominator  of  the  right-hand  side  of  Eq.  (B-  1) 
approaches  the  limit 


n 

E  a  p(u)/9.) 
i=l 

p(ui)  /9_) 


2 

p(u)  /e_)d|a(u)) 


0.~O 

1 

i=l,  .  .  .  ,m 


m 


E 

i,  k=l 


J- 


9  lnp(io  /0_) 

W. 

i 


a  lnp(uu/0_) 


p(u)  /0  )d|i (ud) 


m  n 

+  2  E  E  X.a 
i  =  l  k=m+2  1  k 


a  inp(uu  /e) 

aeT 

1 


7-k)  p(uj  /9_)d|i  (uu) 
p(u)  /0) 
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n  f  pCuo/0.)  p(oj  /e,) 

+  i  - S^eT  P(”/i)d4(“) 


where  —  denotes  partial  differentation  with  respect  to  the  i  component 
of  0.  This1  last  expression  can  be  written  in  a  compact  form  by  introducing 
the  matrices 


D= 


A 

m  x  m 

’"a'” 


m  x  (n-  1-m) 

- IS . 

(n-l-m)  x  (n-l-m) 


A  = 


f  blnp 

J 


lnp(uu  /0_) 


d  lnp(a)/£) 

- g- -  p(u)/0_)d|i(u)) 


A= 


/ 


a  EtaW .  p(»/e)duto) 


ae. 


p(oi/e ) 


i  =  1  ,  .  .  .  ,  m 
k  =  m+1 ,  ....  n-  1 


B  = 


I- 


p^/-i+l) 

p(o)/e ) 


p(^/ik+1) 

-piiTgr 

i,  k=m+l,  .  .  .  ,  n-  1 


and  the  vector 


t  [\ •  •  •  j  \  >  ^ _ i -i i  .  .  .  ,  a  J 

—  1  m  m+2  n 
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With  this  notation,  the  limiting  form  of  the  denominator  of  the  right-hand 
side  of  Eq.  (B-l)  is  D  l.  It  can  be  shown  that  if  g(G0  is  not  constant  and 
if  there  exists  an  unbiased,  finite- variance  estimator  of  g(0_),  then  the 
denominator  of  Eq.  (B-l)  can  never  vanish.  It  follows  from  this  fact  that 
the  matrix  D  is  always  positive  definite. 

With  the  special  parameter  values  chosen,  the  numerator  of  the  right- 
hand  side  of  Eq.  (B-l)  becomes, 


We  have  now  obtained  the  following  form  of  the  Barankin  bound, 


y  'E(0_  -  e)(0  -  0)'  y  * 


a'  ni)2 

V~D T 


(B-2) 


valid  for  all  (n- 1  )-dimensional  vectors _t  and  all  choices  of  6k,  i=m+2,  .  .  .  , 
n.  We  now  maximize  the  right-hand  side  of  Eq.  (B-2)  with  respect  to  L. 
This  can  be  done  via  the  Schwartz  inequality  as  follows, 


(y'  N^)2 
Tm — 


i  i  . 
-  T  ~2  i 
(y'ND  D  l) 

V  D  l 


y_'  ND^N'j 


Therefore, 
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(B-3 ) 


y'  E(0_  -  0J(6  -9)'  j;2y'  N  D^N'y; 


or,  equivalently 


E(9_  -  £)(6  -  9_)'  s  N  D”  1  N/ 
for  all  choices  of  9.,  i=m+2,  .  .  ,  n. 

—l 

The  relationship  between  the  bound  given  by  Eq.  (B-3)  and  the  Cramer- 
Rao  bound  can  be  brought  out  by  means  of  Frobenius1  formula  for  the 
inverse  of  a  partitioned  matrix.  Applied  to  the  D  matrix  this  formula  reads, 


D 


-1 


Noting  that 
where 

we  can  write 


A' 

1  +  A_1A  A_1A' A"1 

l-A-1/ 

-a-'a'a-1 

:  a-1 

i 

'A'1 

N  = 

[mLm  *] 

H 

-  m+2 *  •  '  ’  — n”— 1 

i' 

N  D”  1N/  =  A’1  +  A"  XAa"  1A/A"  1  -  A^Aa'V 
-  $  A~  1A/  A"1  +  §A_1$  ' 

=  A-1  +  (§  -  A_1A)  A-  1  (§  -  A"  *A) ' 


S 
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Therefore,  our  bound  can  be  written  in  the  final  form, 


* 


E(0_  -  e_)(0  -  0_)'  S  A-1  +  ($  .  A_1A)  A”  1  ($  -  A_1A)'  (B-4) 

If  we  pick  n  =  m+1,  the  second  term  on  the  right-hand  side  of  Eq. 
(B-4)  vanishes  and  we  obtain  the  Cramer-Rao  bound. 

E(0  -  0)  (6  -  0)'  ^  A'1 


For  n  >  m+1,  we  obtain  other  bounds  that  are,  in  general,  an  improvement 
on  the  Cramer-Rao  bound.  This  follows  from  the  fact  that  positive  definite¬ 
ness  of  D  implies  positive  definiteness  of  A,  which,  in  turn,  implies  that  the 
matrix  (§  -  A~^A)  A  ($  -  A”  A)'  is  at  least  positive  semi- definite. 


T' 
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APPENDIX  C 


In  this  section  we  will  evaluate  the  Barankin  Bounds,  Equations  (3-3) 
to  (3-8)  for  the  linear  FM  waveform.  We  begin  by  computing  the  A-matrix, 
Eq.  (3-3),  where 


d  s(t;9  ) 

1  _ ° 

Aij  =  N-  -  §9^ 

i,  j,  =  1,  2 


Bs(t;9  ) 
— o 


se . 
j 


dt 


(C-l) 


We  use  the  notation  0Q  to  denote  the  true  value  of  the  parameter  9_.  The 
parameter  vector  is  given  by  9^  =  T,  9^  =  Y  arid  the  signal,  s(t;9_)  is 


s(t;9  )  =  a(t-T)  cosfoi  t  +  i  (t-  t)^  +  y] 

c  ^ 


(C-2) 


so  that 


8s(t;9o) 

§T 


a(t-T  )cos[uu  t  +  ^-(t-T  )2  +  Y  ] 
o  L  c  2  o  o 


Ss(t;  0  ) 
'  -o 

3y 


+  u(t-To)a(t-To)sin[aJct  +  £  (t-TQ)2  +  yJ 
-  a(t-To)sin[u)ct  +  ^(t-To)2  +Yo] 


(C-3) 


(C-4) 


As  an  example  of  the  type  of  manipulations  involved,  we  calculate  the  first 
term  in  detail 


N  A 
o 


11 
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* 


Using  standard  trigonometric  identities  we  get  terms  at  DC  and  at  20)^.  The 
latter  terms  can  be  neglected  and  we  obtain 


NoAll  =  2  I  [a<t-T0))2  dt  I  (t.T 0)2a2(t-T  )dt 

_oo  _  a-> 

(C-6) 

Theoretically  the  modulation  a(t)  is  usually  assumed  to  be  0  for  |t|  >  T/2  and 
V^E/ T  for  |t|  =  T/2  in  which  case  the  first  integral  above  does  not  exist. 
However,  in  practice  such  a  discontinuity  cannot  occur  due  to  the  band  limitation 
imposed  by  the  transmitting  and  receiving  equipment.  In  this  case  a(t)  will  be 
a  "smooth"  function  and  the  contribution  from  the  first  term  would  be  negligible 
compared  to  that  of  the  second.  We  shall  therefore  ignore  its  effects. 
Evaluating  the  second  integral  we  obtain 


A 


11 


2  2 
U  T 

12 


(C-7) 


Using  similar  manipulations  it  is  easy  to  show  that 


12 

"  A21  "  0 

(C-8) 

E 

22 

N 

o 

(C-9) 

Because  of  the  simple  structure  of  the  A-matrix,  it  is  easy  to  calculate  its 
inverse  which  is 


(C-  10) 


Next,  we  compute  the  A-matrix  defined  in  Eq.  (3-4)  as  follows: 
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(C- 11) 


A 


Jj 


J 

_  00 


9s(t;0  ) 

— - 


s(t;  0_Q)]dt 


A0 .  = 

2j  N 
J  o 


r 

_  00 


3 s (t;  0  ) 

5y~°  [s(t;a.)-  s(t;6o)]dt 


(C-12) 


We  use  a_.  to  denote  the  value  of  the  parameter  9  other  than  its  true  value. 
First  we  evaluate  the  term 


r®  3s(t;9  ) 

O 


3t 


s  (t;  0Q)dt 


.r 


{- 


a(t-To)a(t. 


2  r  II 

■  T  )cOS  lu)  t  +  t  (t-T 
O  L  C  Z 


.h 


Vol 


+  u(t-To)a2(t-To)sin[uJct+  j  (t-Tof+Y0]cos[iUct+^(t-Tof+Yo]^dt 

=  0  (C-  13) 

a  result  which  follows  from  the  fact  that  a(t)  is  assumed  to  be  an  even  function. 
Therefore 


N  A 
o 


ij 


ds(t;0  ) 

— 5t  °  s  (t;0Cj  )dt 


Using  the  trigonometric  identities 


B 

(TTTT^Jdt 

(C- 14) 


cos  A  cos  B  =  [cos(A+B)  4-  cos(A-B)] 
sin  A  sin  B  =  Jj-  [sin(A+B)  +  sin(A-B)] 


(C-  15) 
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we  find  that 


A  =  -  1  f  a(t-T.)a(t-T  )cos[ut(T.-T  )-  i  (t  2-  T  2)  +  (y  -Y.)]dt 
o  lj  2  Jf  o'  1  j  o'  2  j  o  v’o  j/J 

00 

p  ^  ^ 

+  7  .  (t-T  )a(t-T  )a(t-T  )sin[a(T  -T  )  -  ^(t  -  T  )  +(Y  -Y.)]dt 

^  CD  O  °  J  J  °  ^  j  °  °J 


(C- 16) 


It  is  instructive  to  evaluate  the  second  integral  first.  Using  the  fact  that  a(t) 
rises  very  quickly  to  the  va lue  V 2 E / T  ,  we  can  use  the  following  approximations 


Case  1: 


T.  -  T  ST 
J  O 


a(t-T.)a(t-r  )  »  0  for  all  t 
j  o 


Case  2:  T.-T  >0,  t.-t  <T 

jo  1  J  o1 


a(t-T . )a(t-T  ) 
J  o' 


2E/T  -  T/2 +  t^  <  t  <  T/2 +Tq 
0  otherwise 


Case  3: 


T.-T  <  0  *  1  T.  -  T  |  <  T 
JO  J  o 


a(t-T.)a(t-To)  « 


2E/T  -  T/2 +tq  <  t  <  T/2 +  t^ 

0  otherwise 


If  we  integrate  term  2,  T^,  between  the  limits  x  and  y,  then  we  obtain 
r-Y 


N 


,t2  =T[lt-V!i*,(Y,ol'ilV’»"lvo'V|dt 

{(t-T0)C0S[ut(T.-T0)-  2(Tj2-T02)  +  (Y0-Yj)] 


X 

-  E 


T(t.-t  ) 
J  ° 


- 7 — - T  sin[ut(T.-T  )-  71  (T  2-T  2)  +  (y  -Y-)]  } 

U(Tj-To)  1  J  o’  2  J  o'  wO  Tj'J  Jx 
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Using  the  appropriate  limits  for  each  of  the  three  cases  we  see  that  for  Case  1, 

N  T~ ,  =  0.  For  Case  2,  we  setx  =  -  T/2  +  t.  and  y  =  T/2  +t  so  that 
o  21  j  y  o 


NoT22  T(T 


coS[u(^  +  to)(t.-to)-|(t.2-t^)  +  (yo-y.)J 


E  (t 

.-t  )  12 

J  o  ( 

-  (--T+T.-T  )  COs[u(-  -^  +  T.)(T.-T  )-  H  (T  2-  T  2)+ (Y  -y.)]> 
2  j  o'  2  J  J  o'  2  '  J  o  'To  'j,JJ 

+  E  7  {si 
uT(t.-t  )  ( 


si„[u(|  +  ToHTj-To)-  |(T.2-T2)  +  (Y0-Yj)] 


sin[u(-  |  +  t.)(t.-to)  -  |(Tj2-To)  +  (V0-Y:i]| 


(C-  17) 


Expanding  the  angles  which  appear  in  this  equation  we  obtain 


U  (-y  +  T  )(T.-T  )  -  H-(T2-T2)  +  (Y  -Y-)  =  (Y  -Y-  )  +  y(T.-T  )  [T-  (T  .  -  T  )] 
v2  o  j  o  2'  j  o  o  j  o  j  2'  j  o  L  j  o 

U(--?  +  T  )(T.-T  )  -  y(T  2-T  2)  +  (Y  -Y-)  =  (Y  -Y-)  -  y(T.-T  )[T-(T.-T  )] 
'2  o  '  j  o  2'  j  o  o  j  o  j  2'  j  o  L  j  o 


We  now  employ  the  definitions 


a  =  (Y0-Yj) 

3  =  £(T.-T  )[T  -  (T.-T  )] 
2  j  o  j  o  J 


and  rewrite  the  preceeding  equations  as  follows 


N  T 


-  E 


o  22  2(t.-T  ) 

J  o 


[cos(a+3)+  cos  (a- 3)]  +  -=  cos  (a- 3)  +- 


uT(t.-t  ) 
J  ° 


2  [sin(a+3 )  -  sin(a  -3 )] 


(C- 18) 


Then  if  we  use  the  trigonometric  identities 
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cos(a+3)  +  cos(a-3)  =  2  cos  a  cos  3 
sin(a+3)  -  sin(a-3)  =  2  cos  a  sin  3 


(C-  19) 


we  obtain 


NoT22  e°s(Vo-Vj)cosl|(Tj-To)[T-(T.-To)] 

+  f  cos{'VYj)  -  J<Tj-To)[T-(Tj-To))} 

+  - — — - 5  cos(y  -y.)sin[ii(T.-T  )[T-(t.-t  )]] 

uT(t  _x  f  °  J  2  J  °  J  ° 

j  o 


(C-20) 


We  repeat  the  same  calculation  for  Case  3  by  setting  x  =  -  T/2  +  tq  and 

y  =  T/2  +  T .  so  that 
1  J 


NoT23  T(t 


)cos[u(^  +  T  )(T  -To).  j(T.2-  T  2)  +tfn-y.)] 


— — - 1  f  (-J  +  T .  -  T 

rvr2  j  ° 

+  T  coslfl-T  +  T^T.-^)  _  |,t.2-To2)+  (Yq" Yj )l| 
UT(T  -T  )  ( 


O  J 


Sin[U(£+T.)(T.-To)-  j(Tj2-T2)+  (YQ-Yj  )] 


-  sin[u(-  j  +  To)(T..To)  -  j(Tj2-To2)  +  (Y0-Yj)]| 


(C-2 1) 


Expanding  the  angles,  we  get 

U(|  +  T.)(T.-To)-  |(T2-T2)+(Y0-Yj)MY0-Yj)+|(Tj-To)[T  +  (T.-To)] 

U(-|  +  T0)(Tj-T0)-  ^(Tj2-T2)+(Y0-Yj)  =  (Y0-Yj)-^TrTo)[T  +  (T.-To)] 


and  now  we  define 
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a  =  (Y0  -  Yj) 

e  =i<yTo>[T+<yTo>] 


and  rewrite  the  preceeding  equations  to  get 


NoT23  =  [cos(a+3)  +  cos(a-P)]  -  ^  cos(a+3) 


J  ° 
E 


uT(t.-t  ) 
J  o 


[sin(a+3)  -  sin(a-3)] 


(C-22) 


Making  use  of  the  identities  used  for  Case  2,  we  obtain 


N  T 


-E 


j.,,  -  ~i - r  cosly  -y.)  cos[^(t.-t  )[T-|t.-t  111 

o  23  (Tj_T0)  °  J  2V  j  o/L  j  oIJJ 


-  f  cos[(yo-Y.)  +  |(t.-To)[T-  |t.-to|]] 


— - Z  cos (Y  -Y. )  sin[|(r  x  )(T-  |t  -t  | )] 

T.-T  )  J  J  J 


uT(t.-t  ) 
J  o' 


(C- 23 ) 


There  remains  the  computation  of  the  first  term,  T^,  for  each  of  the  three 
cases.  Again  we  make  use  of  the  fact  that  a(t)  rises  very  quickly  to  the  value 


V  2E/ T  ,  but  this  time  we  use  the  approximations: 


Case  1: 


T.-T  ^  T 
J  q1 


a(t-T.)a(t-T  )  =  0  for  all  t 
J  o' 
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Case  2 ; 


T.-T  >  0  It.-T  |  <  T 
J  o  ,  J  o' 


•  /  H 

a(t-T  )a(t-TQ)  w  — 6(t-TQ-T/2) 


Case  3:  T.-T  <  0  It.-T  I  <  T 

jo  >  J  o' 


a(t-Tj)a(t-To)  «  -y~  6(t-TQ+  T/2) 

where  6(t)  denotes  the  Dirac  delta  function.  Now  if  we  integrate  term  1, 
Eq.  (C_l6),  for  Case  1  we  get  N  T^j  =  0,  for  Case  2  we  get 


N  T 
o 


coS[u(|  +  To)(Tj-To)-2(Tj2-T02)  +  (V0-Yj)] 


=  +fcOsKY0-Yj)  +  ^(T.-To)[T-(T.-To)]) 
and  for  Case  3  we  get 

NoT13  =-fcosfU(-l  +  To»(Tj-To,-2(Tj2-To)  +  <Yo-Yj)1 
=  -  fcoS{(Yo-Yj)  -  7(T.-T0)[T-  |t.-To|]} 


(C-24) 


(C-25) 


Finally  we  combine  terms  1  and  2  which  for  Case  1  yield  A 
we  combine  Eqs.  (C-20)  and  (C_24)  to  get 


=  0. 


For  Case  2 


E/N 


^  cos(y  -Y.)  cos (T.-T  )[T-  (T  -T  )]] 


(Tj-To)  O-'y—'Z-'  J  O-  J  O' 

E/N 

+  -T-£tcOs[(Yo-Yj)+|(Tj-To)[T-(Tj-To)]] 


+  cos[(Yo-Yj)  -|(T.-To)[T-(T.-To)]]} 

2E/N 

+  — - - -2  cos(Yo-Yj)sin{l(T  -T  )(T-  (T  -T  )]} 

LIT  T.-T  )  J  J  J 

J  o 
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Using  the  identities  in  Eq.  (C-19),  we  obtain  the  final  result  for  Case  2, 


T.-T  >  0, 
J  ° 


ij  no|(t.--, 


cos(y  -Y.)cos{y  (T  -t  )[T  -  (t.-t  )]} 


_  r  i  t  -  t  .  /  (  T  ' 

-T  J  °  J  2  ■  j  o  ‘  '  J  o' 


+  -^,cos(Yo-Yj)  cos[|  (T^-To)[T  -  (Tj-TQ)]} 


+ 


cos(y  -Y-)  sin{y  (T  -T  )[l  -  (T 


uT(t.-t  ) 
J  o 


o  J 


J  o 


rfTo'1}| 


(C-26) 


In  a  similar  way  for  Case  3,  Tj_T0  <  0,  we  combine  Eqs.  (C-23)  and  (C-25)  to 
get 


lj 


•  =  -Jr-  <  /  ~  1 — r  cos(Y  -Y.)cos{y(T.-T  )[T-  |t.-T  I]} 

3  No  (  (Tj-To)  °  J  2  J  °  J  ° 


-  -=■  cos(y  -Y-)cos{y  (t.-t  )[T  -  1  t.-t  I  ]} 


°  J 


J  °  JO’ 


+- 


cos(Y  -Y.)  sin{y  (T.-T  )[T  -  I  T.-T  J  ]}! 


UT(t.-t  ) 
J  ° 


2  o  j'  L2  '  j  o'1  j  o 


(C-27) 


The  second  row  of  the  A-matrix  is  defined  by 


pco  ds(t;  9  ) 

i=TT  J  - [s(t;a.)  -  s(t;9J]  dt 


L2j  N  J  9y 
J  o  -°° 


-J 


— o 


As  before  we  first  evaluate  the  term 


r 

-00 


9  s  (t;  9  ) 


— o 


9y 


f®  ^  ry  » 

a  (t-To)sin[u)ct+2-(t-To)  +  Y0]cos[u)ct+  j(t- Tq)  +YQ]dt 


=  0 
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since  we  are  ignoring  the  effects  of  the  terms  at  2(1)^.  Therefore 


r00  9s(t;0  ) 

N  A,.  =  I  — - — 

o  2j  9y 


r 

s(t;a.)dt  =  -  a(t-T^)a(t-T.)cos[a)ct+2-  (t-T  ^  +  Y„] 


“J 


J 

_  00 


sin[u)  t  +  ^-(t-T.)  +  y.]  dt 

c  2  J  J 


and  using  the  identities  in  Eq.  (C-15),  this  reduces  to 


N 


1  P  ®  ^ 

A  =-T  |  a(t-T  )a(t-T.)sin[ut(T.-T  )-  ^-(t  ,-t  )  +  (y  -y.)]dt 

)  2j  2  J  o'  j'  J  o  2V  j  o  o  ]/J 


for  I  T.-T  I  >  T 

*i  o' 

*  y 

UTIT^TJ  «s[ut(T  -To).  ^(T.2-To2)  +  (V„-Yj )] 


0 

+E 


J  o 


x 


We  perform  the  evaluation  for  the  case  where  Tj_T0  >  0  by  setting  x  =  -T/2+t\, 
y  =  T/2  +  Tq  so  that 


N  A_  .  =  + 


o  2j  uT(t.-Tq) 


[cos[u(^  +  Tj(T.-Tj-  ^(T.2-T^)  +  (Y^-Y,)] 


o'  J  o  L  j  o  O  ] 


-  costul-  l+Tj*Tj-To>  -  7<Tj Z-To>  +  (VYj)1! 

UT(tE-T  )  sin(Y0-Yj>  sin[7(Tj-TolT  -  «Tj-To»] 


J  ° 


(C-  28) 


a  result  which  follows  from  the  manipulations  performed  for  A^.  When 

T.-T  <  0  we  set  x  =  -  T/2  +  T  ,  y  =  T/2  +  T.  and  obtain 
jo  o  j 


NoA2j  =  hT(t.ETo)  sin(Y0-Yj )Sin§  (T . -T0>[T+(t . -to>]) 


(C-2  9) 
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Equations  (C-28)  and  (C-29)  can  be  combined  into  the  single  expression  for 

It.-t  I  <  T 
1  J  o' 


•  2E 


uT(t.-t  ) 
J  O 


sin(Y  -Y  )sin{y(T  -T  ) [ T-  |  T.-T  |  ]} 


(C-30) 


The  B-matrix  was  defined  in  Eq.  (3-5)  as  follows 


Bi  •  =  exp  ^  I  [s (tjOCj^)  -  s(t;£o)][s(t;a.)  -  s(t;9^)]dt 
J  o  -00  ^ 


(C-  3 1) 


First  we  evaluate  the  term 


r 


j 

_  00 


2  2 

s(t;a.)  s(t;a.)dt  =  1  a(t-T.)a(t-T.)cos[u)  t +y(t-T.) +Y.]cos[uu  t+^(t-T.)  +Y.]dt 
— 1_ J  cc  1  J  c  c  i  1  c  ^  J  J 


00 

=  T  a(t-T  )a(t-T  )cos[ut(T  -T  )  -  y(T  2-T  Z)  +  (Y.-V.)]dt 
^  JC0  1  J  J  1  ^  J  1  1  J 


u, -2  _  2, 


UXIT.-T.)  ■inMtTj-Ti)-  Z(T.  -T.  )  t  Iv.-Y. 


)] 


X 


We  make  use  of  manipulations  which  have  already  been  performed  to  write 
the  expression  as 


j  sftiCuMtia.Jdt  =^T"(f~_T  )  cos(Yi-Y-)sin{^(T.-T.)(T  -  |t-T.|)} 
_  oo  J  *  j  o  J  J  J 


The  other  terms  in  the  exponent  can  be  evaluated  from  this  quantity  to  give 
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B  .  .  -  exp  -r-p-  <  1  -  =  . - r-  cus  i  v  •  -  Y  l  s  ini  t  , 

1J  No  (  UTO^-T  )  1  °  1  o 


cos(>VYjsin[y  (t.-tJ(T  -  |t.-to|)] 


—^rr - .  cos (y •  ~ Y  )sin[^(T.-T  )(T-|t.-T  |  )] 

uT(t.-t  )  v’j  o'  1 2  '  j  oM  1  j  ol/J 


J  o 


+  uT(t:-t~)  cos(Yj-Yi)  sin[^(T.-T.)(T-  1  T.-t,  |  )]] 


J  1 


J  1 


(C-3Z) 


The  last  quantity  needed  in  the  evaluation  is  the  $ -matrix  and  this  is 
given  simply  by 


(C-33) 


Some  simplification  can  be  achieved  by  manipulating  the  expression  which 
keeps  recurring  in  all  of  the  formulae,  namely 


7<yTo>[T-  lTj-Tol]  =  Hl--V^  ['  -  |-V^|] 


For  the  linear  FM  signal,  u  =  2-rrAf/T  where  Af  represents  the  signal  band- 

2 

width  and  therefore  uT  =  2-rrAf  T  where  T  Af  is  the  time  bandwidth  product. 

Let  us  make  the  definition 


3  =  2irAf.  T 

so  that  the  key  angle  can  be  written  as 


Similarly  we  can  write 


(C-  34) 


(C-35) 
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and 


(C-  36) 


|iT 

~r 


(T.-T 

J  ° 


1 


T.-T 

J  o 


± 


2 

T 


_2 

T 


1 

+  2(t.-T  ) 
J  ° 

T 


and  finally 


UT 

T 


( 


T.-T 

J 


(C-  37) 


* 


We  next  define  the  normalized  delay  difference 


At. 


J 


(C-  38) 


and  phase  difference 


Ay.  =  Y. -Y 
J  J  ° 


(C-39) 


and  use  the  above  relations  to  rewrite  the  matrices  in  the  following  more 
compact  notation: 


A.. 


J_  _E_ 
T  N 


a. . 
!J 


i  1 , 2,  j  1,2 . n 


(C-40a) 


where 


li  =  (±2  +  a?”)  cosAY.  •  cos[|  At  (1  -  |AT  |)]  +  B — — ^  cosAy.  sin[|  At  (1  -  |  At.  |)] 
J  \  j  /  J  J  J  -^-(AT.)  j  J  J 


1  8 

i2-  =  -  g -  sinAy.  sin[£  At.(1  -  1  AT.  |  )] 

J  -5-at.  j  j 


(C-40bl 


(C-  40c ) 


where  the  +  sign  corresponds  to  At^  >  0. 
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n 


(C-41a) 


B.. 


where 


b..  =  1  -  g1  cosAvi  ‘  sin[^AT.(l  -  1  AtJ  )]  -  -yp —  cos  Ay.  sin[|-AT.(l-  |  At.  |)] 

J  ■£■  At  .  -S-At.  j  J  J 

^  1  2  1 


+  -5 - - -  cos(Ay.-AY.)  sin[y(AT.-AT  )(1  -  j  AT .  -  At.  |  )] 

7  (At. -At.)  1  J  z  1  J  1  J 

2  i  J 


It  was  stated  in  Eq.  (3-6)  that  the  error  covariance  matrix  was  lower 
bounded  by 


E  =  A-1  +  (*  -  A_1A)A_1(§  -  A_1A)  # 


A  =  B  -  A  '  A-1  A 


(C-42) 


- 1  2  2 
The  A  -matrix  is  given  by  Eq.(C-lO),  but  since  \l  T 

rewrite  that  equation  to  give 


=  -^  (ut 

T 


2 


2 

) 


we  can 


A 


-I 


1 

E/N 

o 


(C-43 ) 


Then 
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Then 


12 

7 

12 

7 


a.,a. 0+a_ , a. 


a,.a,„  +a„a. 


12 


-a,  a,-+a.  a. 


12 

7 

12 

7 


12 

7 


and  therefore 


A’A"lA  =  N0(‘2-l-l'+-2-2') 

th 

where  a,  ,  a0  are  vectors  whose  i  rows  are  the  elements  a,  a 
— 1  — Z  li 

Therefore  the  A-matrix  is  given  by 

A  =  B  ‘  W  +  ±2±z) 

o  x  3  7 

an  expression  which  depends  only  on  the  parameters  E/Nq  and  3, 
Next  we  consider  the  matrix 

(*  -  A_1A)..  -  i=l,2;  j  =  l,2,  .  .  .  ,  n 


11  In  Z 1  Z  n 
l  _ a  +a^^a^ 

12  Id  22  Zn 

a.  a.  +a~  a~ 
In  In  Z  n  Z 

2^  respectively. 

(C-45) 
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From  Eqs.  (C-33)  and  (C-44)  we  see  that 


z,  ,  1 2  T 

*lj  =  (Tj-To>  •  alj 

^  =  (Yj-Yo»  -  a2j 


(C-46) 


It  was  stated  that  the  lower  bound  on  the  error  covariance  matrix  was 
given  by 


E(9  -  9  )(  9  -  9  )'  =  E 
—  — o  —  — o 


This  means  that  for  every  vector  t, 


i  = -T  E  - 


In  particular,  when  9  =  (t,  y)/  and  when  l  =  (1,  0)  ,  we  get  a  lower  bound  on 

the  mean- squared  delay  error  when  the  channel  is  incoherent.  Applied  to 
Eq.  (C-42)  using  Eqs.  (C-43),  (C-45),  and  (C-46),  we  find  that 


n  n 


_  / *  x2  .  1  T  ,  ^  V  c 

E(t-TJ  i  pTST  L  L  I’li  6ij  'Pij 


E/N 


o  3 


i=l  i=l 


(C-47) 


th  1 

where  6..  represents  the  ij  element  of  A  .  We  can  normalize  this  result 

1J  2 

with  respect  to  the  time  scale  by  dividing  by  T  ,  which  leads  to  the  expression 


N  N 


1  ,2  ,,  1  12  y  V  ^li  ,  '‘’l 

•tz+  L  L  —  6ii  — 


tzE(?"t°)  e 


i=l  j  =  l 


where  now 


Let  us  set 


T.-T 
1  o 


T  ^  1  i  T 


12 

? 


71  an 


Acpi  =4cPH  -  At.  -  au 

P 


(C- 48) 
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so  that  we  obtain  the  following  lower  bound  for  the  normalized  delay  error  0  ^  : 


_  2  >  1 
CTt  =  E/N 

o 


12 

7 


n 


+  ) 


n 


y 


i= 1  i=l 


Acp.  6..  Acp. 
1  J 


(C 


where  Acfh  is  given  by  Eq.  (C-48),  =  (A  where  A  is  given  by  Eq. 

The  latter  matrix  involves  a,  av.  b..  which  are  obtained  from  Eqs. 

li  2i  ij 

and  (C-41)  respectively. 


-49) 

( C-45). 
(C-40) 
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